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A modified Chapman function with a variable scale height 
gradient has been found to be in good agreement with the electron 
density distribution obtained experimentally within the height 
range of about 100 km below the F2 peak to an altitude of about 
700 km. The scale height distribution derived from this model 
is also consistent with the neutral gas scale height of the 1961 
COSPAR International Reference Atmosphere. 
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ELECTRON DENSITY DISTRIBUTION 
in THF UPPER F-REGION* 


by 

S. Chandra^ 

Goddard Space Flight Center 


INTRODUCTION 

Although ionospheric investigations have been carried on 

the region above the F peak and to obtain some reliable observational data. 

On the basis of their ion trap experiment from Wallops Island, Virginia, Hanson and McKibbin 
(Reference 1) concluded that the ion-density profile between 350 a* 750 km 

static distribution with a constant temperature of 1240°K obt ai„ed 

o * 4fi PM about two hours after sunset. Jackson and Bauer ^Keiere ; 
ment corresponded to 8.40 PM, a f the same sta tion which supported 

rS- Again assuming an ion mass o, 15, they 

concluded the daytime temperature to be about 1640 K. 

The uurpose of this work is to examine electron density data in the light of the above findings and 

L?r"r^ 1,16 observationai data and 

permits the determination of the scale height and other relevant parameters. 


EMPIRICAL MODEL OF THE F-REGION 


Wright (Reference 3) proposed a model of the ionosphere above the ^ 

a simple Chapman function with a constant scaie height o, 100 km “J ~1 *, the 
ence 4), who developed an extensive theory of F2 region, showed that such a simpliiied 


♦Published in substantially the same form in /. 6 eopbys. Res. 68, April 
fNAS-NASA Post-Doctoral Associate. 
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ionosphere was not consistent with experimental results in general. They attempted to fit an experi- 
mental profile obtained by Garriott (Reference 5) with a constant scale height gradient and found close 
agreement only up to 100 km above the F2 peak; the calculated and observed values were found to 
differ considerably in the higher region. Yonezawa and Takahashi suggested that a better fit with the 
experimental results might be obtained by assuming a variable scale height gradient. 

With the availability of more reliable experimental data, it is now possible to examine these 
queshons in more detail. In arriving at a model of the electron density profiles, from near the peak 

o e ayer to the region above, the following known features of the electron density distribution 
have been taken into account: 


. The region above 250 km is essentially isothermal and may be divided into three parts ac- 
cording to the relative concentrations of the ionic constituents. It is now generally believed 
that oxygen ions predominate up to an altitude of about 900 km, helium ions from 900 km to 
about 3000 km, and protons above 3000 km. The thicknesses of these regions vary with 
atmospheric temperature (Reference 6), but the various ionic constituents are in thermo- 
dynamic equilibrium with the neutral atmosphere. 
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The electron distribution near the peak of the F layer can be approximated by a parabolic 
layer and that above the peak by an exponential layer varying as exp -z/h' , where H' is the 
scale height of electron ion gas. This region is usually called an isothermal diffusive equi- 
librium region. Since the predominant ionic constituent in this region is o + and the neutral 

atmosphere also consists predominantly of atomic oxygen, we may write, assuming thermo- 
dynamic equilibrium, 


H' = 2H 0 


( 1 ) 


where H 0 is the scale height corresponding to atomic oxygen. 

In the region near the F peak the scale height of the neutral atmosphere is height dependent 
because, even though the region is isothermal, the mean molecular mass changes with height. 
It approaches, h 0 , the scale height of atomic oxygen at heights well above the peak. 

We will begin with a Chapman function modified to include a variable scale height: 


N m exp 


f h dz l f h dz' 

‘ 1 -H- - H TT 


( 2 ) 


where n corresponds to electron density at height h and h. is the height of the maximum electron 
density, N„. The physical significance of the scale height h, will be discussed later. The following 
analytical form for H has been found to be compatible with the aforementioned features of the F region- 
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where 



h - h 


H „ 


a 


H 0 -H m 

H 0 


and H m is the "scale height" at the peak (z = 0) . Equation 3 has been so formulated that with the 
increase of z' , H approaches H 0 and reduces to a constant scale height for a = 0 . For az '«l , Equa- 
tion 3 assumes a linear form. Substituting Equation 3 in Equation 1 gives 


1 Z ' 

t 

z 

exp TV - ■ i 

1 “ a exp 1 

(-T-) 

~ exp 

i 

1 - a exp 1 

(-T-) 


( 4 ) 


It is easy to verify that Equation 4 reduces to the appropriate forms near the peak and well above the 
peak. In view of the different ion transition regions, Equations 3 and 4 are assumed to be valid a few 
hundred kilometers above the peak where o + ions are predominant. 


COMPARISON OF EXPERIMENTAL AND EMPIRICAL DATA 

From Equation 4, N/N m may be plotted as a function of z', for arbitrary values of a. This is shown 
in Figure 1 for values of a = 0, 0.3, 0.4, 0.5 and 0.6. The curve corresponding to a = 0 is equivalent 
to the simple Chapman function. The parameter a may be interpreted as an index which measures 
the departure from the simple Chapman function. 

In order to verify that Equation 4 is in agreement with the measurements, we shall consider the 
electron density profile obtained by Jackson and Bauer (Reference 2) and the ion-density profile by 
Hanson and McKibbin (Reference 1), shown in Figures 2 and 3 respectively. Equation 3 is completely 
determined for a given set of values (N m , h m ,a and h 0 ) which can be evaluated for any given experi- 
mental profile. The numerical values corresponding to the profiles shown are given in Table 1. The 
computation was performed on the Goddard Space Flight Center IBM 7090 computer' by a differential 
correction least- square fit. 

In Figures 2 and 3 Equation 4 is compared with the experimental results. It is seen that the 
empirical formulation of electron density distribution is in excellent agreement with the experimental 
data. Note that the value of a is higher in the day than in the evening. Since a is an index which 
measures the departure of Equation 4 from the Chapman function, it may be concluded that the electron 
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ELECTRON DENSITY (cm ’ 3 ) 

Figure 2— Comparison of an Experimental Electron Density 
Distribution (Reference 2) with the Empirical Model. 




• HANSON & McKIBBlN (Nov 9 t 1960 ; 2044 EST) 
EMPIRICAL MODEL (a - .39; Ho - 77 km) 
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density thstri*ition tends towards the ~ 

reached by Yonexawa, Martya, Duncan and Dungey (References 

cal considerations. considered three 

— X'U sun conditions, to test tbe ya- 

lidity of Equation 4. 

J +v»o Wallni 


Table 1— Numerical Values 
of the Parameters of Figures 2 and 3. 


Profile 


Electron Density* 
Ion Density t 


♦Reference 2 
tReference 1 


Time 

(EST) 

N 

m 

(1 0 5 cm' 3 ) 

h m 

(km) 

(km) 

a 

1500 

2040 

5.55 

3.73 

280 

313 

102 

77 

0.48 

0.39 


In all three cases the Wallops 
Island Data (taken at 0947 and 2143 
EST, July 13, I960 and 1817 EST, 
June 24, 1961) were in excellent 
agreement with Equation 4 between 
the height range of about 100 km 
below the peak to an altitude of about 

700 km. 


SCALE HEIGHT DISTRIBUTION 


Alter the values ol h. , . and „ 0 are obtained lor a = olile, Z' 

pute H as a lunction ol height. The height va e heights ol the COSPAR International Refer- 

Bauer, Hanson and McKibbin profiles an en atomic oxygen in the region 

ence Atmosphere (Reference 13). Since H approaches the scale g 
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800 



Figure 4-Comparison of Neutral Scale Heights of Cl RA 
(wild curves) with the Scale Heights Derived from Elec- 
tron-Density Data (dashed curves*). 


well above the peak, the corresponding isothermal 
temperature can be computed from the following 
relation 


H - — 

mg ■ ( 5 ) 

where k is the Boltzmann constant, m the mean 
molecular mass, T the temperature and g the 
acceleration due to gravity. By taking m = 16 
and the appropriate value of g at 600 km, the iso- 
thermal temperatures corresponding to the 
Jackson and Bauer and Hanson and McKibbin 
profiles are about 1620 °K and 1180 °K, respec- 
tively. These values are within the limits as- 
signed by these authors to their temperature 
estimates. For the sake of proper comparison, 
the numerical values of the neutral scale heights 
of the reference atmosphere have been normal- 
ized to the temperatures in the isothermal 
region derived from the experimental data. 


K is seen from Figure 4, that the agreement between the neutral scale heights of CIRA and the 

Tr'r r r ron densi,y data is - — — « - 

tron de rT theref0re be reasonable to assume that the scale height derived from the elec- 

ron density profile is actually the scale height of the neutral atmosphere. In the subsequent 
cussion we shall use h for the neutral scale height. suosequent dis- 


EQUATION OF CONTINUITY 


unTeT“! ned ln iy ‘ 8 “ maS ' 

the influence of the gravitional field only, the continuity equation may be written: 


dN d / \ 

dt = q - L - 37 (NW d ) 


(6) 


Jo? a r d ^s^^^^ and l0SS ’ res P ec tively , and w D is the upward ambi- 

ron ion gas, given by the following expression (Reference 14): 


-D 


1 

dN 

(\ 

dT 1\ 

N 

dz + 

(t 

dz + 


6 


(7) 




where h' is the scale height of the electron-ion gas,T is the temperature and the diffusion coefficient 


D 


4.5 x 10~ 17 
' n 


T 2 sin 2 I , 


( 8 ) 


where I is the magnetic dip angle and n is the density of the neutral atmosphere. For an isothermal 
region it can be easily shown that: 


n - n exp 


-It) 


0 ) 


where n m is the number density at the height of the maximum. Equation 6, together with Equations 1, 
2,7,8 and 9, may be written in the following form: 


§ ■ [*■-(£ -*♦•)• ft 


(io) 


where 


X = exp^-J-H-; . 

The last term in Equation 10 represents the effects of diffusion and may be computed for a given elec- 
tron density profile and a given value of n B . 

Assuming that the loss term is an attachment type with a loss coefficient /3 given by Ratcliffe et 
al. (Reference 15): 


( fh-300\ 

/3 = 10 4 expl-J so J sec 

between the height interval of 250-350 km, we can estimate the relative importance of the loss and 
diffusion terms. Further, for a quasi-equilibrium condition, it is also possible to estimate the rate 
of production. Thus, by using the values of N m ,h m , a and T corresponding to Figure 2 and assuming 


Table 2 

Numerical Estimates of the Production, Loss, and Diffusion Terms for Wallops Island. 


Altitude 

L(=/3N) 

div NW d 

q = L + div NW d 

(km) 

(cm -3 sec' 1 ) 

(cm -3 sec' 1 ) 

(dm* 3 sec' 1 ) 

255 

124 

18 

142 

280 

84 

37 

121 

306 

47 

36 

83 

331 

25 

29 

54 

357 

14 

22 

36 
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n m - 2 x lOVcm 3 at the peak (h m - 280 km), the numerical estimates of the production, loss, and dif- 
fusion terms can be tabulated for Wallops Island, Virginia (i = 70°). In Table 2, the last column, 
which represents the production, is computed under the assumption of quasi-equilibrium. 

The numerical values of q shown in the last column are in reasonable agreement with the esti- 
mates of the production rate of 0 + given by Watanabe and Hinteregger (Reference 16). 


CONCLUSION 


An analytical form of the electron density distribution has been found which can be made to fit 
observational data, with high accuracy, within the height range of about 100 km below the F2 peak to 
an altitude of about 700 km. This technique may provide a useful tool for matching - in the region 
around the F2 maximum - the bottomside observations with the data from the Topside Sounder Satel- 
lite. In principle it is possible to use Equation 4 to extrapolate the bottomside data of the quiet iono- 
sphere up to an altitude of about 700 km. 

The scale height distribution derived from the electron density data is consistent with the neutral 
scale heights of the COSPAR International Reference Atmosphere. By assuming that Equation 4 is a 
solution of the equation of continuity and using appropriate values for diffusion and loss coefficients, 
the estimate of electron production rate is found to be in reasonable agreement with the values ob- 
tained from solar flux data by Watanabe and Hinteregger (Reference 16). Consequently, it may be 
possible by further refinement of this method and definitive rocket flights containing simultaneous 
ionospheric, neutral atmosphere and solar radiation experiments, to use ionospheric data for studying 
the neutral atmosphere. 
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